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Abstract 

The gluon fragmentation function of g — > h c at leading order of strong coupling constant a s and 
typical velocity v is calculated in the framework of NRQCD, of which the contributions from both 
color-singlet and -octet processes are taken into account and hence the result is infrared safe. It 
is found that the fragmentation probability of a high energy virtual gluon to split into h c is about 
7.1 x 10~ 7 . The analytic fragmentation function presented in this work can be employed to the 
study of h c production in large transverse momentum at hadron colliders, such as Large Hadron 
Collider (LHC), and is also applicable for the study of hb hadronproduction. 
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It is well-known that the fragmentation mechanism dominates the heavy quarkonium 
hadroproduction at large transverse momentum l|, |2] . It is hence important to obtain the 
corresponding fragmentation function in order to properly estimate the production rate of a 
specific charmonium state. Fortunately it was found that these fragmentation functions for 
heavy quarkonium production are analytically calculable by virtue of perturbative quantum 
chromodynamics(QCD), with limited universal (phenomenological) parameters. 

In past years, Braaten et al. calculated the fragmentation functions of gluon to the S- 
wave charmonium states rj c and J / ib . Later, the gluon fragmentation to polarized J/ib 
was obtained by fitting 141, and the fragmentation function of gluon to J lib via color-octet 
m ech anism wa S ca M a t ld to nex t - t o-lea di n g (N LO ) order ™c y g The tr a„ 10 n 
functions of gluon to P-wave spin-triplet charmonium states Xcj were found infrared cut 
dependent {(j], and after including the color-octet mechanism they then become infrared safe 
Q]. The polarized fragmentation functions of gluon to polarized \ c j were also obtained \]\. 
The fragmentation functions for gluon to D-wave spin-singlet and -triplet were presented 



in [8|, [9]. The fragmentation functions of heavy quark splitting into S-wave charmonium 



states were calculated in 



3. 
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12j, and the corresponding calculations for higher excited 



states, the P-wave and D-wave states, were given in Refs.[13l. Il4j. 1 151. Il6l|. Notice that since 
in this letter we are not going to give a complete review of the calculation of quarkonium 
fragmentation functions, readers who are interested in the details on this respect should 



refer to recent review articles, like [1 

From above brief introduction, one may notice that the fragmentation function of gluon 
splitting into the spin-singlet P-wave state h c is still absent. Because the production rate 
for h c is low and the decay modes of h c are obscure, which hinder experimenters to measure 
it, for a long time the study of h c hadroproduction in large transverse momentum is not so 
urgent and hence the fragmentation function for it. However, in recently the h c is observed 
in experiment [l8, lj], 20, 2l|, and next people may further study its detailed natures. 
With the development of relevant experiments, people expect that more data on h c will be 



collected, like in hadron collision, the LHC for example, in the near future [22j, [23|, [24]. In 
this sense, the fragmentation function of gluon splitting into h c is necessary and valuable. In 
this paper we calculate it at the leading order of a s and v in the framework of non-relativistic 
QCD(NRQCD) 25j, |26(, including contributions from both color-singlet and color-octet. 
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Generally speaking, heavy quarkonium hadropro duct ion via fragmentation mechanism 
involves both short- distance and long-distance effects. The differential cross section can be 



expressed convoluted form 



da H (p) ~ y2 / dz(&i(p/z, fi)Di^ H (z, fi) 

„■ JO 



(1) 



This is a factorized formula, in which the short- and long-distance effects have been sep- 
arated, da is the cross section for parton i production, and is perturbative QCD(pQCD) 
calculable due to the energetic parton and large momentum transfer. D^h{z, p) is a process- 
independent fragmentation function, which describes how parton i with invariant mass p 
splitting into hadron H with light-cone momentum fraction z. For light hadrons, like 
and 7i, the fragmentation processes involve low energy scales, such as Aqcd and m^ n , the 
fragmentation function is therefore a non-perturbative quantity and cannot be computed 
within the scope of pQCD. However, for the fragmentation process of a parton splitting into 
heavy quarkonium, though there will be some non-perturbative effects due to hadronization, 
the production of a heavy quark is a short-distance process and hence is pQCD accessible. 
It turns out therefore that the fragmentation function for heavy quarkonium production can 
be further factorized [3}. 

In the framework of non-relativistic QCD, the fragmentation function of a virtual gluon 
splitting into heavy quarkonium H reads 

D a ^ H (z^) = J2dn(z^)(0\O^\0). (2) 

n 

Here, n is the color-spin-orbital quantum number of the heavy quark pair with null relative 
momentum. The short-distance coefficient d n (z,p), describing the production of heavy 
quark pair with appropriate quantum number n, is pQCD calculable. are local four- 
fermion operators in NRQCD, and their vacuum expectation values are proportional to the 
probabilities of heavy quark pairs with quantum number n hadronizing into quarkonium 
states. According to NRQCD, in (j2J) the short distance sector d n (z,p) can be computed 
order-by-order in strong coupling « s (2m c ), and the long distance sector, the matrix elements, 
can be expanded in series of the typical relative velocity v of heavy quarks inside heavy 
quarkonium {3, 26]. In principle the fragmentation function is calculable to any order in oc s 
and v as desired; in practice, normally the NLO results are enough for phenomenological aim. 
It should be mentioned that in recently there are discussions about whether the definition 



(1) 



(2) 






(6) 



FIG. 1: Feynmann diagrams of the color-singlet process g* — > cc^-P^] + g + 5 
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FIG. 2: Feynmann diagrams of the color-octet process g* — ► cc[ Sq ] + g 



of fragmentation function beyond NLO in strong coupling expansion is complete or not in 
the framework of NRQCD 28|, |29(], whereas it has no influence on our study in this work. 

For process of gluon splitting into h c , there are two classes of processes at leading order 
of a s and v. One is the color-singlet process as shown in Figf]], in which cc pair is in spin- 
singlet, color-singlet and P-wave state (denoted by 1 -P 1 < ' 1 ' > ); the other is color-octet process as 
shown in FigJ2l in which the quark pair is in spin-singlet, color-octet and S-wave state( 1 S'f ) ). 
The full fragmentation function for g — > h c then composes of color-singlet and color-octet 
terms, like 

D g ^ hc (z^) = d 1 (z } A)(0\O h ^CP 1 )\0} + d s (z)(0\O h ^( 1 S )\0}(A) } (3) 
where the short-distance coefficients di(z,A) and d$(z) can be obtained through the cal- 



culation of Feynmann diagrams Figfj] and FigfS], respectively. The factorization scale A is 
introduced to separate the effect at short distance of order l/m c from the one at long distance 
of the radius of heavy quarkonium l/(m c v). (0|C^ c ( 1 Pi)|0) and (0|Cg c ( 1 )S'o)|0) are matrix 
elements of NRQCD operators, scaling as m^v 5 and m z c v b according to velocity-scaling rules 
in NRQCD. Their dependence on A can be obtained by renormalization group equations. 



To leading order of a s (A) they are [261 ]. 

A^OlOfrfPOlO) = 

aA ( o|^(^ )|o) = ^^ (oi^CPOlo). (4) 

The last equations can be used to sum up large logarithms of m c /A when the factorization 
scale are chosen to be much smaller than m c [2j. i.e., 

(0|^(^ )|0)(m c ) = (0\O h s ^S )\0)(A)+ J° F \n^-(0\Ot( l Pi)\0), (5) 

with QCD beta function (3 = (11N C — 2nf)/6. In following analytical results, we will 
show that ds(z) is A independent at leading order of a s . Thus, the dependence on A in 
the factorization formula ([H comes from (0|£>8 C (%)|0) and difoA), the same as in the 
fragmentation functions of g — > Xcj- 

Leaving the polarization of fragmenting gluon unsummed, we obtain the squared ampli- 
tudes of color-singlet processes, 

+ (terms proportional to or q v ) . (6) 

Here, p, k and q are the momenta of the cc pair, final and initial gluons, and s = q 2 
is the virtuality of fragmenting gluon. In the condition go 3> m c which is required to 
make fragmentation mechanism a reasonable approximation, the fragmenting gluon is almost 
transversely polarized. And under this condition we can drop terms proportional to <x, or q v , 
which are considered as gauge artifact and suppressed in an appropriate axial gauge The 
Lorentz indices in Eq. flB]) imply that there are angle dependence in the squared amplitudes. 
After averaging the directions of momentum p and k, we get the following replacements, 

p^Pu ^ -^\p±\ 2 g^u, (7) 
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k^k v — > ~9pv 2 "_ ^-\p±\ 2 (3cos 2 9 k - 1) + s(y 2 - r)(l - cos 2 fe ) , (8) 
z y r L z 



2 tcosOk 



kflv+PpK -> , = , (9) 



with 



(l + r)/2-f-y + *y 2 / . 2^ / 1rV > 

cos9 k = , \p±\ = s(-r + 2yz - z ) , (10) 

where the dimensionless variables are defined, 



p • q k ■ q 4m 2 p° + p 3 

y = — , t = , r = , z= (11) 

s s s q v + q 6 



Notice that the squared amplitude of Eq.Q is now proportional to metric tensor g^ u , and 
the fragmenting gluon behaves like an on-shell one. 

To avoid the presence of infrared divergence while gluons in final state become soft, 
we set a lower cutoff A on the energy of the final gluons, then the color-singlet coefficient 
d\(z) converges and depends on A. After integrating over the three-body phase space, the 
color-singlet part of fragmentation function in Eq.Q is obtained, 

dUz.A) = s = I dr I dy- -ry^ / ^ h(y,r) 

u ' ; 20736^^ J (r+z , )/2z y (1 - y)% 2 - rf\ (1 - r) 2 



i=0 



{- [-2z 2 + 3z + 2(1 - z) ln(l - z)\ In — + w(z)\, (12) 



162irm^. 

with functions w(z), h(y,r), fi(y,r) and gi(y,r) read 

w{z) = 2(1 - 2)ln 2 (l - z) + [-2 - 3(-2 + z)z + 3(-l + z)\nz] ln(l - z) 

+ i[27r 2 (l - z) + 3z(-7 + 8z) + 3z(-3 + 2z)ln.z] - 2(1 - z)Li 2 (l - z) - (1 - z)Lj 2 («) , (13) 
6 

h(y, r) = (1 + llr - 5r 2 + r 3 )(l - 20r + 6r 2 - 4r 3 + r 4 ) + (2 - 20r + 606r 2 - 280r 3 + 94r 4 - 20r 5 + 2r 6 )?/ 
+ (4 + 468r - 792r 2 + 232r 3 - 44r 4 + 4r 5 )j/ 2 + (8 - 1120 + 560r 2 - 96r 3 + 8r 4 )j/ 3 

+ (-240 + 1072r - 208r 2 + 16r 3 )j/ 4 + (544 - 448r + 32r 2 )j/ 5 + (-448 + 64r)j/ 6 + 128j/ 7 , (14) 
/o(»,r) = r 3 (36 - 15r + 25r 2 - 57r 3 + 59r 4 - 97r 5 + 23r 6 - 7r 7 + r 8 ) 

+ (-164r 3 + 54r 4 + 80r 5 - 34r 6 + 344r 7 - 2r 8 + 12r 9 - 2r 10 )y 
+ (-82r 2 + 274r 3 - 374r 4 - 46r 5 - 746r 6 + 26r 7 - 14r 8 + 2r 9 )y 2 
+ (-4r + 488r 2 + 88r 3 + 652r 4 + 824r 5 - 564r 6 - 108r 7 )y 3 
+ (1 + 33r - 1277r 2 - 637r 3 - 877r 4 + 1947r 5 + 553r 6 + r 7 )y i 
+ (2 - 160r + 1734r 2 + 488r 3 - 2794r 4 - 1032r 5 + 66r 6 )j/ 5 
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h(V,r) 
h{y,r) 



go(y,r) 



9i(V,r) 
92(y,r) 



-336 + 80r - 1136r 2 - 1264r 3 )y 8 + (448 + 736r + 992r 2 )y 9 + (-384 - 448r)i/ 1 

r 2 (6 + 8r - 47r 2 + 33r 3 - 24r 4 + 46r 5 - 7r 6 + r 7 ) 
-16r 2 + 96r 3 + 56r 4 + 2r 5 - 206r 6 - 10r 7 - 2r 8 )y 



16r 6 )y 3 

-7 - 285r - 592r 2 + 320r 3 + 151r 4 - 99r 5 )y 4 + (34 + 336r + 204r 2 + 128r 3 + 242r 4 )j/ 5 

84 - 220r - 340r 2 - 316r 3 )/ + (112 + 232r + 248r 2 )y 7 + (-96 - 112r)y 8 + 32y 9 ] , 
r 4 (-9 + r - 8r 2 - 8r 3 + r 4 - r 5 - 8r 6 ) + (46r 4 + 10r 5 + 71r 6 + 20r 7 + 4r 8 + 66r 9 + 7r 10 )y 



-28 + 356r - 1144r 2 + 2232r 3 + 644r 4 - 396r 5 )/ + (136 - 352r - 352r 2 + 400r 3 + 968r 4 )j/ 7 



128?/ 



-16r - 148^ 



2A8r° 



- 74r* 



- 502r° 



110r 6 + 2r 7 )y 2 



28r 3 - lllr 4 - 94r 5 - 148r 6 



■ 34r' 



- 213r 8 - 20r 9 )y 2 



-148r 3 + 193r 4 + 146r s + 150r 6 + 356r 7 - lllr 8 - 26r 9 )«/ 3 



192r^ 



166r J 



120r 5 - 330r 6 + 695r 7 + 182r 8 )f/ 4 
1604r 6 - 498r 7 + 16r 8 )y 5 



2 - 590r 3 - 662r 4 + 1930r 5 + 614r 6 - 124r 7 )y 6 
1160r 2 + 1152r 3 - 896r 4 - 112r 5 +408r 6 )y 7 

8r - 1152r 2 - 216r 3 - 696r 4 - 7Mr B )y 8 + (-104r + 360r 2 + 864r 3 + 816r 4 )?/ 9 
8 + 216r - 376r 2 - 520r 3 )j/° + (-16 - 48r + W0r 2 )y 11 ] , 
2y92(y,r) , 

r 3 (-3 - 5r + 8r 2 + 3r 4 + 13r 5 + 16r 6 ) - (-10r 3 - 14r 4 + llr 5 + 4r 6 + 84r 7 + 138r 8 + llr 9 )?/ 
-10r 2 + llr 3 + 114r 4 + 60r 5 - 196r 6 - 495r 7 - 76r 8 )y 2 
16r 2 - 183r 3 - 492r 4 + 86r 5 + 940r 6 + 201r 7 - 8r 8 )y 3 
16r + 175r 2 + 796r 3 + 640r 4 - 914r S - 215r 6 + 62r 7 )?/ 4 
-96r - 688r 2 - 1172r 3 + 108r 4 - 76r 5 - 204r 6 )?/ 5 

218r + 922r 2 + 562r 3 + 502r 4 + 372r 5 )y 6 + (200r + 472r 2 + 568r 3 + 408r 4 )?/ 7 
-4 + 36r + 268r 2 + 260r 3 )j/ 8 - (8 - 8r - 80r 2 )?/ 9 . 



(15) 
(16) 



(17) 



(18) 
(19) 



(20) 



Compared to the color-singlet short-distance coefficient, the calculation for the color- 
octet one is much more straightforward and the result does not depend on A, 

d 8 (z) = ^j-[-2z 2 + 3z + 2(1 - z) ln(l - z)} . (21) 

Now we have the final expression for fragmentation function D g _f lc (z, /i) at the leading order 
in strong coupling constant. If the value for factorization scale is chosen to be much smaller 
than m c , we should exploit the renormalization group equation (J3J), or its solution to 
sum up the logarithm of A/m c . In this work, the factorization scale is set to be m c so that 
we can avoid the large logarithms [2]. The value of the scale \i in a s are chosen to be 2m c , 
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which is the square root of the minimum value of the virtuality of the fragmenting gluon. 

Integrating the fragmentation function D g ^h c (z, 2m c ) over the momentum fraction z, we 
get the fragmentation probability, 

JdzD 3 ^2ra c ) = -A-^ i ^ i <Of.) + (0.33)(*)( mc )] . (22) 

To make the total probability positive, the lower bound for the magnitude of color-octet 
matrix element at factorization scale A = m c is set to be 

(Ot)(m c ) > ^^(^ c > ■ (23) 

For numerical evaluation, the input parameters are taken as follows: m c = 
m hc /2 = 1.78GeV, a s (2m c ) = 0.26, (Oie^Pi) |0) = 0.32GeV 5 which enables 



O^SoWirric) > 6.3 x 10- 3 GeV 3 , (0|O^( 1 5 )|0)(m c ) = 9.8 x l(T 3 GeV 3 [24], |3lJ. With 



above inputs the fragmentation probability in Eq. fT22l is found to be about 7.1 x 10 7 . Com- 
pared with the probabilities of gluon fragmenting into y c0 , y c i, and y c2 , which are 0.4 x 10~ 4 , 

fl 

1.8 x 10 4 and 2.4 x 10 4 [2|, the value of gluon fragmentation to ^t c is smaller by two to 
three orders, and it is even less than the color-singlet probability of process g* — > J / tpgg by 
an order [3]. The z dependence of the fragmentation function at /i = 2m c and A = m c is 
shown in Figure [31 

In conclusion, we have computed in this work the fragmentation function of gluon to 
P-wave spin singlet quarkonium h c . The analytic expression is infrared safe while both 
color-singlet and -octet processes are taken into account. It is found that the fragmentation 
probability of a high energy virtual gluon splitting into h c is about 7.1 x 10~ 7 . The obtained 
analytic fragmentation function can be employed to the study of h c hadroproduction, espe- 
cially in the large transverse momentum region. Finally, it is worthy to note that result in 
this work can be readily applied to the study of hj, physics. 
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FIG. 3: The fragmentation function D g ^f lc (z, fi) as a function of z for \i = 2m c and A = m c , the 
dotted line is the contribution from color-singlet process, the dash-dotted line is the contribution 
from color-octet process and the solid line is the total result. 
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